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Abstract 



I \ We present a detailed numerical study of the one-dimensional Holstein model 

c/3 , with a view to understanding the self-trapping process of electrons or exci- 



tons in crystals with short-range particle-lattice interactions. Applying a 



'^ \ very efficient variational Lanczos method, we are able to analyze the ground- 

O \ state properties of the system in the weak- and strong-coupling, adiabatic 



and non-adiabatic regimes on lattices large enough to eliminate finite-size 
effects. In particular, we obtain the complete phase diagram and comment on 



> 

00 

fSj , the existence of a critical length for self-trapping in spatially restricted one- 

(N 



dimensional systems. In order to characterize large and small polaron states 
we calculate self-consistently the lattice distortions and the particle-phonon 
correlation functions. In the strong-coupling case, two distinct types of small 



'^ ' polaron states are shown to be possible according to the relative importance 



of static displacement field and dynamic polaron effects. Special emphasis is 
on the intermediate coupling regime, which we also study by means of direct 



o 
o 

> 

X. 

H diagonalization preserving the full dynamics and quantum nature of phonons. 

The crossover from large to small polarons shows up in a strong decrease of 
the kinetic energy accompanied by a substantial change in the optical absorp- 
tion spectra. We show that our numerical results in all important limiting 
cases reveal an excellent agreement with both analytical perturbation theory 
predictions and very recent density matrix renormalization group data. 



I. INTRODUCTION 

Electrons, holes or excitons delocalized in a perfect rigid lattice can be "trapped" within 
a potential well produced by displacements of atoms from their carrier-free equilibrium 
positions, provided the particle-lattice interaction is sufficiently strong ^^ . Trapping of a 
carrier in this manner is more advantageous energetically as compared to wide-band Bloch 
states, if the lowering of the carrier's energy due to its binding exceeds the strain energy 
required to produce the trap. Since the potential that binds the carrier depends on the 
carrier's state itself, i.e., the local distortion of the lattice is self-induced by the particle, 
this process is called "self-trapping" or "auto-localization". Obviously, self-trapping (ST) is 
an highly non-linear phenomenon. A self-trapped state is referred to as "large" if it extends 
over multiple lattice sites. Alternatively, if the quasiparticle is practically confined to a single 
site the ST state is designated as "small". Nonetheless, ST does not imply a breaking of 
translational invariance. ST eigenstates in a crystal are Bloch-like. Thus coherent transport 
of ST particles can, in principle, occur but the width of the corresponding band is extremely 
small (cf. the discussion in 0). 

Introducing the concept of polarons into physics, the possibility of electron immobi- 
lization or ST was pointed out by Landau as early as 1933 |^. Shortly after the ST of 
excitons was also suggested and studied theoretically 0,|§]. ST polarons and excitons can 
be found in a wide variety of alkali metal halides, alkaline earth halides, II-IV and group 
IV semiconductors, condensed rare gases, organic molecular crystals, electrochromics, and 
other systems P,§. With the observation of polaronic effects in new materials exhibiting 



exceptional properties such as the high-T^ cuprates [10| or the colossal magnetoresistance 



manganates 0], research on polarons has attracted renewed attention. 



Although the problems of exciton and electron ST have much in common there are 
fundamental differences. Most notably excitons are short-living non-equilibrium quasipar- 
ticles being immediately after the optical excitation in the free state and can reach the ST 
state only by tunneling through the potential barrier at low temperatures. Moreover, if 



the electron and hole, forming the exciton, have very different effective masses, the internal 
coordinates of this large-radius (Wannier-Mott) exciton will be of importance. 

It is clear that the microscopic structure of the ST state is very diverse in various groups 
of materials. The stability of different types of ST states depends on the nature of the 
electron/exciton-phonon (EP) coupling (e.g. deformation potential |T2| against Frohlich |TB 



coupling), the vibrational frequencies (e.g. acoustic vs. optic), the dimensionality (D) of the 
lattice and other parameters. A detailed classification of ST states and ST criteria is pre- 
sented, e.g., in the excellent review of Rashba 0. In particular, from a scaling treatment of a 
continuum lattice model in the adiabatic limit, it has been shown that in multi-dimensional 
systems {D > 1) with only short-range non-polar EP interaction a ST carrier always forms 



a small polaron |jTj,|T5|,^. If there is an energy barrier that separates delocalized and spatial 
confined states, the free and ST states can coexist |16| , p!7| , as evidenced, e.g. for solid xenon, 
by the coexistence of two exciton luminescence bands (one narrow and the other broad) . 
On the other hand, the presence of a long-range polar EP coupling ensures that at least the 
formation of large polaron states with moderate lattice deformation becomes possible. The 
picture is qualitatively changed when turning to low-dimensional systems. Recently it was 
shown that, unlike the continuum limit, the formation of a ST state in the 2D case within 
a model of local EP coupling is always accompanied by the formation of an energy barrier 
attributed to the lattice discreteness ||18[. The ID case is essentially different from the 2D 
and 3D ones. In ID ST proceeds without overcoming a barrier. To be more precise, the ST 
state is the ground state of the infinite system at any value of the EP coupling. In the weak 
EP coupling regime its radius exceeds several times the lattice spacing, i.e., a large polaron 
is formed even under the non-polar interaction condition [0]. 

Although the basic concepts underlying the ST transition are long standing and the 
gross features of "large" and "small" polarons and excitons have been extensively stud- 
ied, our understanding of the ST problem is still incomplete. In particular the physically 
most important crossover regime, characterized by intermediate EP coupling strengths and 
phonon frequencies, is difficult to handle theoretically due to the failure of the standard 



phase transition concept [0,^. As yet, there exist no well-controlled analytical techniques 
to describe the transition region. Other problems, for example, concern the existence of a 
critical length for ST in spatially restricted ID systems E^], the behaviour of the polaron 



kinetic energy [^], or the spectrum of light excitons under ST conditions 0. 

With these motivations, in this paper we want to discuss the ST problem using numerical 
methods. The focus is on the Holstein model in one- dimension. By exact diagonalization of 
finite systems we analyze various ground state and spectral properties of the model. Since 
for the Holstein model exact results are very rare ||20| , |23|j2^ and previous numerical studies 



have been limited either to small systems or to a particular parameter regime p5|-|29|l this is 



a challenge by itself. Besides, we hope to gain more insight into the physics of small/large 
polarons and into the nature of the localization-delocalization transition. 

The paper is organized as follows: In the next section we briefly introduce the Holstein 
model and outline our variational Lanczos approach that allows to study the ground state 
properties for all regimes of parameters on large lattices in a very efficient way. The numerical 
results will be presented in Sec. HI. More precisely, the phase diagram of the Holstein model 
(A), the electron lattice correlations (B), the optical response (C) and the kinetic energy 
(D) will be discussed. The principal results are summarized in Sec. IV. 

II. MODEL AND METHODS 

With our focus on the self-trapping phenomenon in systems with only short-range non- 
polar electron- or exciton-lattice interaction, it is appropriate to consider a model, where 
a single excess carrier is placed in an one-dimensional periodic array of identical molecular 
units, each having an internal vibrational degree of freedom. Introducing electron (exciton) 
af and phonon b^ destruction [creation] operators we can write Holsteins's molecular crystal 
model |12[ in lattice site representation as 



i i (ij) 



(n^ = ala^] below ^ = 1). In the case of electrons the Holstein Hamiltonian (1) has been 
studied extensively as a paradigmatic model for small polaron formation. Here t denotes the 
nearest-neighbour free-electron transfer amplitude, Ep is the strong-coupling polaron binding 
energy (in the atomic limit t = 0), and uq is the bare phonon frequency of the dispersionsless 
optical phonon mode. Referring to excitons we have in mind the small-radius (Frenkel or 
charge transfer) excitons only, and neglect, in the lowest order of approximation, the internal 
structure of the exciton, i.e., we consider it as a single neutral particle. 

Generally speaking, the ground-state and spectral properties of the model (1) are gov- 
erned by three ratios (control parameters) defined from the bare energy scales t, ujq, and Ep. 
First the adiabaticity parameter 

a = ujo/t (2) 

determines which of the two subsystems, excitons/electrons or phonons, is the fast and 
which is the slow one. In the adiabatic limit a <^ 1, the motion of the particle is affected 
by quasi-static lattice deformations (adiabatic potential surface) depending on the strength 
of the particle-phonon interaction. On the contrary, in the anti-adiabatic limit a ^ 1, 
the lattice deformation is presumed to adjust instantaneously to the position of the carrier. 
Conveniently the particle is referred to as a "light" or "heavy" electron/exciton in the adia- 
batic or non-adiabatic regimes 0. The second parameter is the dimensionless EP coupling 
constant 

X = ep/W, (3) 

where W = 2Dt denotes the half-width of the electron (exciton) band in a rigid D 
dimensional lattice. Let us stress that A represents the ratio between "localization" energy 
(oc Ep) and bare kinetic energy {W) of a single particle. Both, A and a, are commonly used 
as parameters within a perturbative analysis of the Holstein model in the limits of weak 
(A ^ 1) and strong (A ^ 1) EP couplings. In the latter case two different approaches, 
based on expansions in powers of {a <^ 1) and (1/A ^ 1), have been elaborated for the 



adiabatic [12,16| and non-adiabatic pO] regimes, respectively. A third parameter. 



5-^ = ep/c^o , (4) 

will show to be crucial in the strong-coupling situation, g"^ determines the relative deforma- 
tion of the lattice which surrounds the particle. 

In the limit of small particle density, a crossover between essentially delocalized carriers 
and quasi-localized particles is known to occur from early quantum Monte Carlo calcula- 
tions [^, provided that the two conditions A ~ 1 and (7 ~ 1 are fulfilled. So, while the first 



condition is more restrictive if a is small, i.e. in the adiabatic case, the formation of a small 



ST state will be determined by the second criterion in the anti-adiabatic regime p9| , |32| , |33 

It is not surprising, that the standard perturbative techniques are less able to describe 
the system close to the crossover region, where the energy scales are not well separated 
(A ~ 1, (7 ~ 1). Therefore we will apply in the following two distinct numerical methods 
that allow investigating the ST phenomenon on finite clusters with great accuracy. 

The first method is a variational Lanczos technique (developed originally for the Holstein 
t-J model [^,^]), which enables us to study the ground-state properties of fairly large sys- 
tems. In recent work this technique has been adapted successfully to treat the lattice degrees 
of freedom in the generalized double-exchange Hamiltonian commonly used for the descrip- 



tion of colossal magnetoresistance materials ^5|. In the case of the pure Holstein model, as 



a first step, we perform an inhomogeneous variational Lang-Firsov (IVLF) transformation 

n=U^ni{, W = e-'^i(^') e--^^^^) e-'^^(") (5) 

with 

S,{A,) = -^Y.Mbl-b,), (6) 

52(7) = -igY^i^l - h) 4«i , (7) 

i 

53(r) = -lnr-V2^(6j6l-6A), (8) 

i 

(rescaling Ti = Ti/t and measuring, in what follows, all energies in units of t). In a cer- 
tain sense the canonical transformation W is a (variational) synthesis of two different ap- 
proaches developed in the adiabatic |]6|, |36| , p!2| , p!6| and anti-adiabatic 0J|,0 theories of ST 



polarons/excitons. Si introduces a set of static site-dependent displacement fields Z\j related 
to local lattice distortions. This transformation ensures the correct behaviour in the adia- 
batic limit. That is, Si describes the ST of "light" excitons/electrons under the conditions 
that (i) the electronic bandwidth significantly exceeds the characteristic phonon frequency 
and (ii) the lattice deformation energy is large (which allows one to treat the lattice vibra- 
tions quasi-classically) . Within polaron theory such type of ST quasiparticle is often called 
adiabatic Holstein polaron (AHP) ||12|. S2, on the other hand, describes the ST process in 
the anti-adiabatic limit. The variational parameter 7 (with < 7 < 1) is introduced as 
a measure of the non-adiabatic phonon dressing of "heavy" particles, designated as "local- 



ized" excitons ||3^ or non-adiabatic Lang-Firsov small polarons (NLFP). For 7 = 1, the 



well-known Lang-Firsov displaced-oscillator transformation results. In addition, we have 



applied the two-phonon squeezing transformation ^a (0 < r < 1) |]39| , ^ . The squeezing 
phenomenon is a many particle effect being of special importance at intermediate EP cou- 
pling strengths. This effect can be seen as a phonon frequency softening and tends to offset 
the (polaron) band narrowing. As a second step, we approximate the eigenstates \^) of Ti 
by the variational states 

|<^y) = \%h) ® \^ei) . (9) 

Then, performing in a third step the average over the transformed phonon vacuum, 7Y = 
{\Ppf^\7i\\Ppf^) , we obtain the effective (electronic/excitonic) Hamiltonian 

-2A7(2 - 7) E^. - e-^'"'^' E(4s- + 4«.) • (10) 

In (10), the first term leads to an increase of the zero-point energy of the phonons if r^ < 1. 
The second and third contributions give the elastic energy and the particle-lattice interac- 
tion, respectively, both owing to the static lattice deformation. As a result of the incom- 
plete LF transformation we get a constant (polaronic) level shift and an exponential band 
renormalization (fourth and fifth terms). Even the simplified model (10) cannot be solved 



exactly. Therefore we carry out a Lanczos diagonalization on finite A^-site lattices using 
periodic boundary conditions. Employing the Hellmann-Feynman theorem, the A^ + 2 vari- 
ational parameters can be obtained by iteratively solving the following set of self-consistency 
equations: 

A = 4A(l-7)n., (11) 

r' = , , (13) 

where 



Hi = in, 



)^ with ^ni = l, (14) 



Ekin = - ^((al«i + «i«i))w , (15) 

EAn = ^{Ani)ff_, (16) 

i 

denote the local particle density, the kinetic energy, and the EP interaction contribution 
to the ground state energy, respectively. Note that each iteration step involves the exact 
diagonalization of 7^(7, r^, {Z\i}). We note further that the Hamiltonian (10) potentially 
contains symmetry-broken states which originate from inhomogeneous displacement fields 
Z\j 7^ 0. Therefore we have to work with an unsymmetrized set of basis states. 

The second method we are going to use in the computational work is the direct numerical 
diagonalization of the initial Holstein Hamiltonian (1). On the one hand this should bring 
out valuable information on the applicability of various approximative analytical and nu- 
merical approaches. In particular we would like to test the quality of IVLF-Lanczos scheme 
described so far. On the other hand, combining our exact diagonalization (ED) algorithm 
with the Chebyshev recursion and maximum entropy methods |^, we are able to discuss 
dynamical properties of the systems, e.g. the optical conductivity, in more detail. Moreover, 
ED provides the only reliable tool for treating the transition regime. Differently from the 
IVLF-Lanczos treatment the translational invariance of the system is ensured. The ED 

8 



method, based on a well-controlled truncation procedure of the phonon Hilbert space, has 



been described elsewhere P2| , |32| , |33[] . Using parallel computers, we are able to diagonalize 
systems with a total dimension of 10^°. 

III. NUMERICAL RESULTS AND DISCUSSION 

A. Phase diagram 

In the numerical work we start with a discussion of the ground-state properties of the 
transformed Hamiltonian (10). Applying the IVLF-Lanczos technique presented in the 
previous section, we have determined the phase diagram of the ID Holstein model. The 
results are depicted in Fig. 1. First let us consider the regimes II and III, corresponding 
to large and small polarons, respectively. Just for brevity we will use in the following the 
"polaron terminology", keeping in mind that all statements hold for the case of Frenkel 
excitons as well. The distinct types of polarons, found in II and III, may be characterized 
by the spatial extension and strength of the (inhomogeneous) lattice displacements Z\j and by 
the magnitude of the polaron variational parameter 7 (see below). From exact analytical |2^ 



and numerical [^,^,^ (cf. also Sec. Ill D) results it is well known that the large-size 
polaron turns continuously into a small-size polaron with increasing EP coupling. Since 
there is no true phase transition between large and small ST states at a > 0, the transition 
line (stars) shown in Fig. 1 only indicates the crossover region, which gets wider as the 
phonon frequency (a) increases. Within our IVLF-Lanczos treatment, the transition line 
has been fixed by the criterion Ai/Aq = 1/e. Performing a finite-size analysis of the II^III 
transition, we found that the results obtained for the 64-site lattice almost agree with the 
extrapolated values for the infinite system. Via Eq. (11), Z\j is directly related to the polaron 
density at site i. Of course, this condition can only give a crude estimate of the "transition" 
from large to small polarons, in particular in the non-adiabatic regime where the static Ai 
are less significant. According to the importance of the Ai {a < 1) and 7 (a > 1) effects. 



the small polaron formed in region III will be called adiabatic Holstein polaron (AHP) and 
non-adiabatic Lang-Firsov polarons (NLFP), respectively. 

As a peculiarity of our finite system a further strongly finite-size dependent transition 
to a nearly free electron state (I) is observed by lowering the EP coupling strength. In 
other words, it seems that a critical coupling Xc{N) or equivalently a critical system size 
Nc{X) exists for self-trapping in ID. Indeed, for the ID continuum model, where the ST 
problem can be described by a non-linear Schrodinger equation, it has been shown recently 



by Rashba pll that the ST condition is 

A>Ac = vrV2A^. (17) 



This relation holds rigorously within the adiabatic theory [^ and is reproduced by our lattice 



model calculation as well (cf. inset Fig. 1). At a = 0, the nearly free electron phase (I) 
corresponds to a solution with 7 = and Ai = A = 4X/N. The kinetic energy, however, 
is unrenormalized. Our IVLF-Lanczos scheme allows to extend the above considerations to 
the finite phonon frequency regime. Again, at low EP couplings, we found a nearly free 
electron phase with a small uniform level shift (oc — 2A[7(2 — 7) + 2(1 — 7)^/A^]). More 
significantly, since we have 7 > now, the inclusion of non-adiabatic phonon effects slightly 
renormalizes the electron bandwidth W = 4exp{ — (7^7^r^}. If A becomes larger than Ac(a) 
the ST proceeds by a monotonic lowering of the total energy without overcoming a ST 
barrier. The scaling of Ac (a) with A^ is shown in the inset of Fig. 1. In the thermodynamic 
limit A^ -^ 00 we obtain Ac(a) -^ 0, i.e., in an infinite ID system ST takes place at any 
finite value of the EP coupling. 

To elucidate in more detail the different nature of polaronic states occurring in the 
ground-state phase diagram of the effective model (10), we present in Fig. 2 the behaviour 
of the variational parameters. 

First of all, we should emphasize that, studying the single-electron problem, the squeez- 
ing effect (r^ < 1) is only of minor importance. This is obvious from (13), which yields 
r^ = 1 in the thermodynamic limit. For finite systems, the leading (l/A^)-corrections 

10 



(oc fyf^7^e~^ "^ "^ Ekin/aN) tends to zero in the weak- and strong-coupling adiabatic and 
anti-adiabatic limits. 

The spatial extension of the static lattice deformation (Z\j) is visualized in Figs. 2 (a) 
and (b) for different EP couplings corresponding to the nearly free, large and small polaron 
cases. As discussed above, the Ai are being constant for A < Ac(a;) (phase I). For the 
large-size polaron (phase II), the lattice displacements fits extremely well to the relation 

Z\i = Z\osech2[Ae//-i], (18) 

where lima^o '^e//('^; ct) = ^- It is worth emphasizing, that the functional form (18), which 
has been derived in the framework of an adiabatic continuum theory [|16],^|18|, also describes 
the displacement fields in the non-adiabatic large polaron regime. Obviously, Ag// defines 
a characteristic inverse length scale in the system, i.e., the radius of the large polaron is 
approximately given by ~ (2Ae//)^^. For a = 0.1 and A = 0.25 we got Xeff/^ — 0.935. 
On the other hand, at a = 3 and A = 2.5, the effective coupling becomes strongly reduced: 
Ae///A ~ 0.116. In the strong-coupling regime, we observe an exponential decay, Z\j ~ 
Z\oe~*'^, of the lattice distortion away from the polaron site (see insets), where ^ denotes the 
small polaron radius. We found ^ ^ 0.29 (0.19) for A = 1.5 (5.0) and a = 0.1 (3.0), i.e., in 
both cases the ST state is mainly confined to a single lattice site. While, in the framework 
of our interpolating theory, the A^ can be taken as a measure of the "adiabatic character" 
of the polaronic quasiparticle, its "non-adiabatic part" is described by the Lang-Firsov 
variational parameter 7 shown in Fig. 2 (c). Of course, in the case of "light" electrons 
{a < 1), the non-adiabatic polaron effect is rather small. In particular for A > 1, when 
the small AHP is formed, 7 becomes strongly suppressed. Here the renormalization of the 
polaron band is mainly driven by the static displacement fields Ai. Otherwise, for "heavy" 
electrons, we observe larger values of 7, which increase with increasing EP coupling strength. 
As a result the free-electron band is transformed into a renormalized polaron band. Due 
to the (generalized) Franck-Condon factor e~^ '^ ^ , the bandwidth is exponentially small 
under strong coupling conditions {g"^, A ^ 1), where also the K-dependent corrections to 

11 



the band dispersion become negligible |3^. Within our variational approach we found a 
first-order transition to the AHP state at extremely strong EP interaction {Xn/m — 4, for 
a = 3). However, it should be noted that this sharp transition is in some sense an artifact 
of our IVLF-scheme that compares the ground-state energies of the AHP and NLFP states, 
both obtained in the lowest order of approximation (remind that by deriving (10) we have 
performed the average over the zero-phonon state only). Including higher-order corrections, 
the NLFP with 7 ^ 1 is stabilized in the non-adiabatic strong-coupling regime (cf. the 
discussion in Sec. HI D). 

B. Electron lattice correlations 

In looking for a characterization of the different polaronic regimes for the quantum- 
phonon Holstein model (1), we have calculated the (normalized) correlation function between 
the electron position i = and the oscillator displacement at site j, 

^'' 2g{no) ' ^''^ 

by means of direct diagonalization technique. Here fermion and boson degrees of freedom 
are related by the well-known relation {b] + bj)'n = 2g{ni)-H (= '^g/N for the single electron 
case) . 

Alternatively, working with the effective Hamiltonian (10), the electron-lattice correla- 
tion function (19) can be expressed as 

Xo,, = l5o, + ^ . (20) 

Hence we can use the static correlation functions xo,j ^i-nd xq,j to test the accuracy of the 
IVLF-scheme. 

Figure 3 shows the static correlation functions (19) and (20) in the adiabatic (a) and 
non-adiabatic (b) regimes for several coupling parameters A corresponding to the different 
polaronic "phases" indicated in Fig. 1. For parameters close to the adiabatic weak-coupling 
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regime (phase I), the amphtude of xo,j is extremely small and the spatial extent of the 
electron-induced lattice deformation is spread over the whole lattice. In the quantum- 
phonon model (1) the coupling gives rise to a weak dressing of the electron at any finite a. 



However, the carrier is not trapped due to the zero-point quantum lattice fluctuations pi 
From (20) it is clear that in the effective model (10) the on-site dynamical polaron and 
spatially extended static displacement field contributions are well separated. Neglecting the 
residual polaron-phonon interaction, the IVLF-approach describes the real situation by a 
(variational) superposition of both effects. In the adiabatic large polaron region II a much 
better description of the exact behaviour is obtained. Especially when the small AHP state 
evolves at A ~ 1 the IVLF results are even in quantitative agreement with the DMRG 
(density matrix renormalization group) data obtained very recently by Jeckelmann and 
White IQ (see inset Fig. 3 (a)). At this point we would like to emphasize that performing 
such DMRG- and, in particular, ED-calculations requires much more memory and CPU- 
time resources than our extreme simple and very fast IVLF-computations. In Fig. 3, the 
system sizes treated within the IVLF-scheme have been restricted in order to make possible 
a direct comparison with the available ED/DMRG data. Obviously, for intermediate to 
strong EP couplings, the IVLF-results agree almost exactly with the ED and DMRG data 
(Fig. 3 (b)). Here the electron-lattice correlations become very localized and finite-size 
effects are less important. Although the behaviour of Xoj? shown in the main part of Fig. 3, is 
found to be very similar for A = 1.5, a = 3.0 and A = 4.5, a = 1.0, we would like to emphasize 
that both parameter sets describe completely different physical situations. The distinct 
nature of the corresponding small polaron states becomes apparent from the variation of 
the static displacement fields shown in the insets. For A = 4.5 and an intermediate (or low) 
phonon frequency, we observe a static lattice distortion in the vicinity of the electron only. 
Since Aq/4X ~ 1, in (20) the second term dominates the first one and we obtain a small AHP 
confined to a single site. Contrary, for A = 1.5 and a = 3, we are still in the large polaron 
region II due to the high phonon frequency (cf. Fig. 1 and the spatial extension of the 
static lattice distortion shown in the left inset of Fig. 3 (b)). Nevertheless the correlations 
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between the electron and the phonon remain local. But now, since the Aj/A\ are small for 
all j, the peak structure of Xoj at j = results from the first term in (20). That means 
it is mainly triggered by the 7-effect (7 ^ 0.69). This interpretation is substantiated by 
our ED results yielding a mean phonon number in the ground state of about 0.625, which 
clearly indicates that the zero-phonon state is still the most probable one. Therefore the 
approximation we applied by deriving (10) and (20) is justified. 

C. Optical response 

Extremely valuable informations on the low-energy excitations in interacting 
electron/exciton-phonon systems can be obtained by studying their optical response. Actu- 
ally, the optical absorption of small polarons is distinguished from that of large (or quasi- 
free) polarons by the shape and the temperature dependence of the absorption bands which 
arise from exciting the ST carrier from or within the potential well that binds it [Q. Fur- 
thermore, as was the case with the ground state properties, the optical spectra of light and 
heavy electrons/excitons differ essentially as well 0. In the most simple weak-coupling 
and non-adiabatic strong-coupling limits, the absorption associated with photoionization 
of Holstein polarons is well understood and the optical conductivity can be analyzed ana- 
lytically p5H49[] (for a detailed discussion of small polaron transport phenomena we refer 
to the review of Firsov 15^]). Serious problems, caused, for example, by the complicated 
behaviour of the adiabatic potential surface, arise if one tries to calculate the spectrum of 
self-trapped light excitons . Moreover, the intermediate coupling and frequency regime is 
as yet practically inaccessible for a rigorous analysis. On the other hand, previous numeri- 
cal studies of the optical absorption in the Holstein model were limited to very small 2- or 
4-site clusters [^,^. In order to calculate the optical conductivity numerically in a wide 
parameter range on fairly large systems, we have implemented our computer code, which 
is based on a combination of Lanczos diagonalization, Chebyshev recursion and maximum 



entropy methods ETi, on parallel machines. 
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The real part of the optical conductivity, Rea{uj) = VS{uj) +a^'^^{uj), can be decomposed 
into the Drude weight (5-function at a; = and a regular part {u > 0) written in spectral 



representation at T = as 51 



a'^Huj] 



\{%\iEj{a]aj^^ - a]_^^aj)\^J\ 



2 



E ' '^ J,' ' 5{uj-E^ + Eo). (21) 



Em — Eq 



In (21), a'^'^^{uj) is given in units of vre^ and we have omitted the 1/A^ prefactor. For the 
discussion of the optical properties it is useful to consider the spectral weight function 



S'^Huj] 



r du'a''''{u') (22) 

"^ U 



as well. 

Numerical results for both, cr^'^^(co') and iS''^^(u;), are presented in Fig. 4. We will start 
with the somewhat more subtle case of light electrons. 

But first let us recall that, restricting ourselves to phononic states \s)ph = 
riilo^ ~Y^ \ i) ' \^)ph '^ith at most M phonons, a ii'-symmetrized state of the Holstein 
model is given as |!P'^) = J2m=o^s=i '^k''^ l-^j"^;^); where S'(m) = {N — I + my. / {N — l)\m\ 
(for more details see P3|). K denotes the toto/ momentum of the coupled EP system. Then, 



if the EP coupling is finite, the ground state \\Pq) and all excited states \^m) contain com- 
ponents that correspond to m-phonon states (with m = Ylf=o^ iT-i ^ M, nf G [0,r7i]) in the 
tensorial product Hilbert space of electronic and phononic states. When the EP coupling 
is small (A ^ 1), these multi-phonon states have less spectral weight, i.e., the phonon dis- 
tribution of the ground state, |c[J^p(M) = J2s\ck=o\'^ , exhibits a pronounced maximum at 
the zero-phonon state [|33|. The maximum of |c™p(M) is shifted to larger values of ?7i as A 



mcreases. 

Keeping this in mind and notice further that in (21) an optical transition can take 
place only within the K = sector (|!^o) — \^k=o))^ ^^e peak structure of o"^^^ shown in 
Fig 4 (a) may be easily understood in connection with the single particle spectrum. For 
low phonon frequencies {a <^ W), the energy to excite one phonon lies inside the bare 
tight-binding band E}^ = —2tcosK and we observe a flattening of the coherent band 
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structure Ek at large momenta [|32| , p2| . Then the coherent bandwidth AE = E.^ — Eq 
is approximately given by a (-C [Ej^ — Eq ]), i.e. by the phonon frequency, where the 
states at finite momenta are predominantly "phononic" states with less "electronic" spectral 
weight. Thus, although in principle an optical excitation can be achieved by "adding" 
phonons with opposite momentum to these states (in order to reach the K = sector), the 
overlap to the mainly "electronic" ground state is extremely small. Therefore we found, 
roughly speaking, the first transitions with non-negligible weight to the free electron states 
and its vibrational satellites (see S^'^^{uj)). This is perfectly illustrated by the inset of 
Fig. 4 (a). Here the first and second group of peaks is approximately located at the bare 
tight-binding energies, E\^ (+ n- a), for the allowed wave vectors of the 8-site lattice used 
in the numerical calculation {K = 7r/4, 7r/2, . . .). Of course, with increasing the lattice size 
the K values will become dense and we will obtain the monotonous decay of the optical 
absorption coefficient observed for large polarons above the photoionization threshold [Q. 
To understand the changes in the optical absorption in the crossover region from a large-size 
polaron to small AHP, the main part of Fig. 4 (a) shows a'"^^ at two intermediate EP coupling 
parameters. In this case the coherent band structure Ek gets stronger renormalized, but, 
more important, the phonon distribution function in the ground state, |c[{^p(M), becomes 
considerably broadened. For instance, at a = 0.1, we have |c™p(M = 25) ^ 0.008 {m = 0), 
0.095 (m = 7), 0.008 (m = 18) and |CP(M = 25) ~ 0.0002 (m = 0), 0.1 (m = 12), 
0.0005 {m = 24) for A = 0.9 and A = 1.0, respectively. Therefore the overlap with excited 
multi-phonon states is enlarged and the optical response is enhanced. The lineshape of 
the absorption bands refiects the phonon distribution in the ground state, where the small 
oscillations are due to the discreteness of the phonon frequency. As a result, the peak 
structure is smeared out and the wide side-bands, belonging to different electronic momenta 
(e.g., K = IT /A and K = tt/2), merge with each other. 

We now turn to the heavy electron case (Fig. 4 (b)). The inset again illustrates the be- 
haviour in the large polaron regime (A = 1, a = 3; cf. Fig. 1), where a rather moderate band 



renormalization occurs (AE ~ 2.33 due to the fiattening effect |^). But now the phonon 
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frequency is large compared to the finite-size gaps between the first mainly "electronic" 
excitations [Ex ~ Eq < a for K < 7r/2). Therefore, in contrast to the light electron case (cf. 
Fig. 4 (a)), the different absorption bands (each built up by several "electronic" i^- levels) 
can be classified according to the number of phonons involved in the optical transition. As 
can be seen from the main part of Fig. 4 (b), the absorption spectrum for a small-size po- 
laron is quite different from that of a large polaron. According to the results of Sec. Ill A, 
a small NLFP is formed in the strong-coupling non-adiabatic limit. Here the phonons will 
heavily dress the electron and, concomitantly, the "electronic character" of the resulting 
strongly renormalized small polaron band becomes suppressed (cf. the discussion of the 
i^-dependent wave-function renormalization factor Z^' = \{^kWk\^o)\'^ ^^ -^^f- El)- -^^^ 
our parameters (A = 6, g"^ = 4), the maximum in phonon distribution function is located 
between m = 3 and 4. The renormalized bandwidth is small compared to all other energy 
scales {AE ~ 0.0782 ^ a, W). Since the current operator connects only states having sub- 
stantial overlap as far as the phononic state is concerned, multi-phonon absorptions (i.e., 
non-diagonal transitions ^B[) become now increasingly important in the optical response. 
This leads to the peak structure observed for the non-adiabatic small polaron optical con- 
ductivity in Fig. 4 (b). Obviously, the different bands are being separated by multiples of 
the bare phonon frequency. The height of the jumps in the u;-integrated conductivity is 
directly related to the probability of the corresponding absorption process. We found that 
substantial spectral weight stays in the lower energy part of the spectrum at frequencies 
comparable to Ep = 2A (^ rrimax ■ «). These absorptions, resembling to some extend a large 
polaron's absorption, are signatures of a ST polaron with intermediate size. In the extreme 
strong-coupling limit the dominant absorption process results from the transfer of the ST 
carrier to the neighbouring site without changing the lattice distortion. That means, the 
optical absorption spectrum should exhibit a single-peak structure at u = 2ep = 4A, which 
corresponds to the lowering of the electronic energy associated with the small-polaron for- 
mation p3| , ^ . This feature already evolves for the coupling strength considered in Fig. 4 (b) 
(cf. AS^^^ for the 7^^ absorption band). 
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D. Kinetic energy 

Further information about the transition from large to small polarons can be obtained 
from the behaviour of the polaron kinetic energy -Efcm- Replacing Tihy 7i, the kinetic energy 
can be easily obtained from static correlation function (15). On the other hand, according 
to the f-sum rule, E^in is directly related to the cu-integrated optical conductivity, 

5*°* = -^ + 5'^"^ . (23) 



2 ~ 27re2 



Calculating, via (21) and (22), iS*"^^ = 5'"^^(oo) numerically, allows us to determine the Drude 
weight T> as well. Sometimes one defines an effective polaronic transfer amplitude p8| , |22 



^p,e// = EkinW / Ekin{0) , in order to characterize the polaron mobility. In our reduced 
units we have tp^eff = '5*°*. From (23) it is obvious that tp^eff includes both coherent and 
incoherent transport processes. Hence tp,e// substantially differs from the exponential factor 

2 

e~^ , obtained in lowest order perturbation theory, and cannot be used to determine the 
coherent bandwidth AE under strong coupling conditions (A, g^ ^ 1). 

The evolution of the kinetic energy (5*°*) as a function of the EP coupling A is displayed 
in Fig. 5 (a) and (b) for the case of light and heavy electrons, respectively. In agreement 



with previous numerical results PT|J28|J2^P^ , E^m clearly shows the crossover from a large 
polaron, characterized by a iS*°* that is only weakly reduced from its non-interacting value 
(5*°* (A = 0) = 1), to a less mobile small AHP/NLFP in the adiabatic/non-adiabatic strong- 
coupling limit. 

For low phonon frequencies {a = 0.1; Fig. 5 (a)), we found a rather narrow transition 
region. As recently pointed out by Capone et al. [^], the decrease of iS*°* in the crossover 



regime is driven by the sharp fall of the Drude weight. By contrast the optical absorption due 
to inelastic scattering processes, described by the regular part of the optical conductivity, 
becomes strongly enhanced ||3^ (see the behaviour of S^'^^). It is worth emphasizing that 
the IVLF-results for 5*°* (A) are in excellent agreement with the ED and DMRG data. 

The large to small polaron transition is considerably broadened at high phonon frequen- 
cies (a = 3.0; Fig. 5 (b)). Here the IVLF-results start to deviate from the exact ones 



when g^ gets much larger than one, thus making the lowest-order zero-phonon approxima- 
tion inherent in the IVLF-scheme less justified. As mentioned already in Sec. Ill A, the 
non- analytic jump-like behaviour at A ~ 4 is an apparent shortcoming of the variational 
approach which compares ground-state energies only. 

Although for large enough g and A the simple formula W = W exp{—g'^}, which should 
not be identified with the "Lang-Firsov approach" p2| , |55[| , works perfectly well in the deter- 



mination of the coherent bandwidth {AE ^ W) |2^j3^, the need of going beyond the lowest 
order of approximation to obtain reliable results for the kinetic energy has been emphasized 
many times |ll|, |33|j4^ , |5^j55[| . In the non-adiabatic strong-coupling limit {g'^ ^ 1, a > 1), 



the ground-state energy obtained within second-order perturbation theory is a tiny little 
bit lower than the IVLF-energy and almost coincides with the ED result. Adapting the 
second-order strong-coupling approach presented in our previous work |3^ to the ID case 
(with 7 = 1), the kinetic energy is obtained, via E^^^"^ = tdt{H), as 



pSCPT__^/h _ -g^ 

^kin — \ / o 2 ^ 



4 / 1 , 

Q, \ 5 / K,=g-' 



(24) 



Here (. . .)« means the average with respect to the Poisson distribution with parameter k. 
As can be seen from Fig. 5, at large EP interactions, the strong-coupling perturbation 
theory (SCPT) gives a sufficiently accurate description of 5*°* in both the adiabatic and 
non-adiabatic regimes. 

IV. SUMMARY 

The main objective of this study was to re-examine in detail the self-trapping transition 
of electrons and excitons in one dimension (ID) within the framework of the Holstein model 
by the use of computational techniques. The calculations are performed by exact diago- 
nalizations of finite systems, where the full dynamics and quantum nature of phonons was 
taken into account. Therefore the results are unbiased and allow to test the applicability 
of a numerically much more efficient variational (IVLF) Lanczos approach proposed by the 
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authors. This IVLF-Lanczos technique is designed to analyze the ground-state properties 
of strongly correlated electron-phonon (EP) models on fairly large lattices, including static 
displacement field, non-adiabatic polaron and squeezing effects. 

Let us summarize the main outcome of our work. 

Fig. 6 illustrates the basic physics contained in the single-particle Holstein model. This 
model describes an continuous transition from large to small polarons with increasing EP 
coupling strength. Depending on the adiabaticity of the system, a = uJo/t, the crossover 
regime is determined by the more stringent of the two conditions A = ep/2Dt and g^ = Sp/uoQ. 
Thus starting from "light" {a < 1) or "heavy" (a > 1) electrons it is possible to understand 
the formation of small adiabatic "Holstein" (AHP) or non-adiabatic "Lang-Firsov" (NLFP) 
polarons as two limiting cases of a general picture. 

In the infinite ID Holstein model the self-trapped state is the ground-state for any 
value of the EP coupling. By contrast, in a finite ID Holstein system a critical length 
or equivalently a critical EP coupling strength exists for self-trapping, rigorously at least 
at cjq = 0. This may be of importance for spatially restricted systems like Cqo- The 
large-size polaron is characterized by spatially extended lattice deformations. Within our 
IVLF-description, the variation of the displacement fields follows the (adiabatic) formula 
Ai oc sech [Ae// ■ i] even in the non-adiabatic regime, but with a strongly reduced inverse 
length scale given by an effective coupling constant Ae//(A,Q;) <^ A. The mean phonon 
number in the large polaron ground state is rather small. In 2D and 3D Holstein systems a 
large polaron phase does not exist. 

The small polaron state is basically a multi-phonon state characterized by strong on-site 
electron-phonon correlations. Due to a large local static lattice distortion the AHP becomes 
quasi-localized on a single site. Also the NLFP is immobile because it has to drag with it 
a large number of phonons in its phonon cloud. 

The transition from large to small polarons is accompanied by significant changes in 
their optical response. Especially for "light" electrons the spectral weight of the regular 
part of the optical conductivity is strongly enhanced at the transition. For "heavy" electrons 
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the lineshape of the optical absorption spectra is highly asymmetric in the intermediate- 
to-strong coupling regime and therefore differs considerably from the usual small polaron 
hopping behaviour obtained for g^ ^ 1. 

As a result of self-trapping the mobility of the charge carriers is reduced. The kinetic 
energy indicates that the crossover region from large to small ST states is rather narrow 
(broad) in the adiabatic (non-adiabatic) regimes. The formation of small adiabatic Holstein 
polarons is accompanied by a dramatic kinetic energy loss mainly driven by a sharp drop 
of the Drude weight. Since both the bandwidth and the (electronic) spectral weight of the 
small polaron band are exponentially reduced in the extreme strong coupling case g"^, A ^ 1, 
coherent small polaron transport becomes rapidly destroyed by thermal fluctuations. 

Finally, we would like to stress that the simple IVLF-Lanczos approach provides an ex- 
ceptional good description of the ground-state properties of the Holstein model, in particular 
for the light electron case. 
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FIGURES 

FIG. 1. Phase diagram of the ID Holstein model. Nearly free polarons, large polarons, 
and small polarons exist in the regions I, II, and III, respectively. Results are obtained for 
finite rings with A^ = 32 (o), 64 (□, *), 96 (A), and 128 (y) sites using the IVLF-Lanczos 
method. The inset shows the critical coupling for self-trapping, Ac, at a = (x), 0.1 {O), 
1.0 (>), and 3.0 (<). The solid line gives the relation (17). For further explanation see text. 

FIG. 2. Variation of the displacement fields Ai away from the central site for several 
characteristic EP couplings A in the adiabatic (a) and non-adiabatic (b) regimes. The 
dependence of the Lang-Firsov polaron variational parameter 7 on A is depicted in (c). 

FIG. 3. Electron-lattice correlations Xo.j in the adiabatic weak-to-intermediate EP cou- 
pling (a) and non-adiabatic intermediate-to-strong EP coupling (b) cases. IVLF-Lanczos 



results are compared with ED data and the DMRG results taken from Ref. [^j . 

FIG. 4. Optical absorption in the ID Holstein model. For N = 8 and M = 25, the 
regular part of the conductivity a^"^^ (thin lines) and the integrated spectral weight S^'^^ 
(thick lines) are plotted as a function of u at different EP couplings for the "light" (a) and 
"heavy" (b) electron cases. 

FIG. 5. Kinetic energy [in units of (— W^)], 5**°*, and contribution of o"™^ to the f-sum 
rule, S^'^^, as a function of the EP coupling. A, in the adiabatic (a) and non-adiabatic 
(b) regimes. Circles (stars) denote exact (DMRG [^) data obtained for a lattice with 
A^ = 8 (32) sites. IVLF-Lanczos results (solid curves) are compared with the predictions 
of standard first-order {W/W = exp{—g'^}; thin solid line) and second order (chain-dashed 
curves) perturbation theory. 

FIG. 6 Schematic phase diagram of the Holstein model. 
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